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I Abstract 

We investigate an antiferromagnetic s=l/2 quantum spin system with anisotropic spin ex- 
change on a fractal lattice, the Sierpihski gasket. We introduce a novel approximative numeri- 
cal method, the configuration selective diagonalization (CSD) and apply this method to a the 
Sierpihski gasket with N=42. Using this and other methods we calculate ground state energies, 
spin gap, spin-spin correlations and specific heat data and conclude that the s=l/2 quantum 
antiferromagnet on the Sierpihski gasket shows a disordered magnetic ground state with a very 
short correlation length of ^ ~ 1 and an, albeit very small, spin gap. This conclusion holds for 
Heisenberg as well a for XY exchange. 
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1 Introduction 



Low-dimensional quantum antiferromagnets (AFM) have been intensively investigated since the de- 
velopment of quantum mechanics in the early twenties PP. A renewed interest was motivated by 
the discovery of high-temperature superconductivity [2j and the peculiar interplay of the magnetic 
and electronic properties of these systems, where antiferromagnetism and superconductivity appear 
in close vicinity (see e.g. and ref. therein). Since then various experimental findings for materi- 
als for electronically one- or two dimensional magnetic systems, such as CuGeO CaVO, SrCuBO, 
have called for more detailed theoretical investigations of the ground state and the low temperature 
properties of one- and two-dimensional quantum AFM (see e.g. j3] and ref. therein). Despite many 
theoretical efforts, many properties of low-dimensional quantum AFM, in particular the interplay of 
quantum fluctuations and magnetic order near quantum critical points, need further explanation. 

The radical difference in the behavior of one- and two-dimensional AFM has been subject of 
current debate, in particular with regard to its interaction with changes in the lattice structure (in 
the cuprates), the presence of spin-peierls transitions or the influence of disorder. One of the most 
significant dimension- dependent properties is the type of magnetic order: the ground state of the one- 
dimensional s=l/2 Heisenberg chain remains disordered P], but for the two-dimensional Heisenberg 
quantum AFM on square, triangular or honeycomb lattices one observes a Neel-like magnetic long 
range order in ground state. The dimensional crossover between d=l and d=2 has been studied via 
investigations of ladder structures or by varying exchange parameters on two-dimensional lattices 
(e.g. on the square lattice or on the triangular lattice) El |Hl El UHl HH II2j- 

In addition to the dimensionality the spin anisotropy can also influence the magnetic order in 
quantum AFM. For example, it is known that in the one-dimensional linear chain an infinitesimal 
small Ising-like exchange anisotropy induces a Neel like magnetic order in the ground state. In zig-zag 
ladders the effects of XY anisotropy may lead to spiral ordering [13j. For two-dimensional lattices 
at finite temperatures (where the Mermin- Wagner theorem forbids any Neel-like long-range order 
for pure Heisenberg exchange) an XY exchange anisotropy can induce a vortex type ordering at the 
Kosterlitz-Thouless transition [Tij . 

We have previously studied the influence of dimensionality on the magnetic order by considering 
a quantum AFM on a particular lattice geometry, the Sierpihski gasket, with a topological dimension 
between one and two. We considered a Heisenberg interaction between the nearest-neighbor spins on 
this lattice and investigated its properties with exact diagonalization and variational wave functions 
[T3t I16j. We supplemented this analysis with thermodynamical properties using the quantum deci- 
mation technique (QDT) [T7] and a decoupled-cell Monte-Carlo approach (DCM) [TH]. Recently we 
extended this research to higher spins [19J as well as to anisotropic spin interactions [201 • From all 
calculated data we have presented arguments in favor of a disordered ground state of the Sierpihski 
gasket quantum AFM. 

All mentioned numerical approaches have some limitations and disadvantages. The exact diag- 
onalization is subject to the well known constraint on the system size. Its application to quantum 
spin models is always challenging but in the case of a fractal lattice even more complicated due to 
loss of translational symmetries. Therefore only small systems with N=15 could be investigated so 
far (In ['2W\ we include a system with N=28, which has a similar topology like N=15 but is in fact not 
a true Sierpihski gasket). The variational wave function calculations suffer from the uncertainties in 
the reference wavefunction due to lattice frustration and the thermodynamical properties calculated 
with QDT and DCM might not probe the true ground state. Therefore all previous conclusions have 
been drawn with particular care and indeed, especially the investigation of larger systems is very 
desirable for further support and verification. 

In this paper we will apply a new technique, the configuration selective diagonalization (CSD), 



2 



in a particular efficient implementation to investigate larger finite lattices. This approach is based 
on ideas developed in quantum chemistry j2Il 1221 evolves around an extrapolative calculation 
of the ground state and low excitations. Using this approach we are able to access larger systems 
in a numerically controlled approximation. We will present results for the N=42 Sierpihski gasket 
(the next larger Sierpihski gasket after N=15) for both the Heisenberg and the XY model. This 
is to our knowledge the largest quantum spin system treated so far with a direct diagonalization 
technique. The treatment of this cluster permits the analysis of spin-spin correlations for larger 
lattice separations and allows us to further strengthen the predictions of a disordered ground state 
for the Heisenberg model. It also enables us to draw similar conclusions for the XY model on the 
Sierpihski gasket. 

In section |2] we will introduce the model and its basic properties. In section El we present all 
numerical methods used in our investigation. We will emphasize on the CSD in very detail, as this 
is the ffist application of this approach to quantum spin systems. In section 0] the results of the 
calculations will be presented and conclusions about the magnetic order behavior will be made. We 
summarize the paper in section O 

2 The Model 

We consider the quantum s=l/2 AFM with anisotropic spin exchange: 

H = JY. (S^SJ + SfSf + AS^S^^). (1) 

<id> 

The antiferromagnetic spin exchange J > is taken between nearest neighbors on the Sierpihski 
gasket (an example of this lattice with N=15 is given in FigC]). The anisotropy A will be studied for 
the Heisenberg model at A = 1 and for the XY model at A = 0. The most important geometrical 
property of the Sierpihski gasket is its fractal Hausdorff dimension oi dj = ~ 1.58. The number 

of spins on this lattice is given by N= |(3" + 3) with n = 1, 2, 3, . . .. In the paper we will focus on 
N=6, 15 and 42 (i.e. n=2,3,4). 

The classical ground state is a planar spin state with 3 sublattices for A < 1. The spins in 
such a sublattice are ferromagnetically aligned, between spins belonging to different sublattices we 
observe an angle of 120°. This ground state is depicted by the arrows in the right part of FigO The 
classical ground state of the Sierpihski gasket is analogous to the ground state of the two-dimensional 
triangular lattice and has no non-trivial degeneracy. 

We would like to mention here that there is an ongoing debate of the relation between the 
properties of the classical ground state and the magnetic order of the quantum ground state. For the 
kagomelaXiice it was argued in [221 ^^^^ infinite degeneracy of the classical ground state is closely 
connected to the high number of low-lying singlets in the quantum case (and therefore maybe to a 
liquid-like quantum ground state.) Ongoing research did show that there are systems like the planar 
pyrochlore lattice with a classical non-trivial degeneracy of the ground state, yet with valence bond 
type long range order in the quantum regime [25. The multi-spin exchange model is another counter 
example with no non-trivial classical degeneracy but with a large number of low lying singlets [211 • 
As it has been discussed already in ^21 the Sierpihski gasket is yet another example as it has no 
non-trivial classical degeneracy, even though in the quantum case it also shows a large number of 
low-lying singlets and a disordered hquid-like ground state. 
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Figure 1: Left: The classical ground state configuration of the N=15 Sierpihski gasket with 3 sub- 
lattices: A (circles), B (squares) and C (triangle). Right: The spin direction in the classical ground 
state with an angle of 120° between spins belonging to different sublattices. 



In this paper we will use a variety of methods to investigate the Sierpihski gasket. For finite quantum 
AFM of small size the exact diagonalization (ED) with a Lanczos approach is the tool of choice to 
investigate the ground state and the low energy spectrum. Because of the exponentially growing 
Hilbert space the approach is usually limited to systems with up to N=36, only for highly symmetric 
lattices like the square lattice one can reach N=40 177] . As already stated the CSD will be used 
for the calculation of the ground state and first excited state of the N=42 Sierpihski gasket. For 
thermodynamic properties we will deploy a complete diagonalization (CD) for smaller finite systems 
and a quantum decimation technique (QDT). 

Whereas the complete diagonalization is routinely used (for more information see e.g. [28]). 
the QDT is not known to a wide audience even though it has been successfully implemented to 
investigate low-temperature thermodynamics of different low-dimensional AFM 129, 30|. Therefore 
we will shortly describe the two basic steps of this approach on the example of the Sierpihski gasket. 
In the first step, one splits the Hamiltonian H of an infinite system into Hamiltonians Hi of finite 
6-spin subsystems (see right part of FigH} {H = J2iHi). Recalling that the renormalization group 
(RG) procedure should preserve the partition function and symmetry of the system and, additionally, 
the decimation procedure should preserve the correlation function, one traces out some spin degrees 
of freedom in each finite subsystem. In the second step one puts these finite subsystems together 
and obtains a renormalized Hamiltonian H': 



Note that while splitting the Hamiltonian into finite subsystems and subsequently replacing the true 
local Hamiltonian by the renormalized H- one neglects the non-commutativity of the spin operators. 
This two-step RG transformation enables one to calculate the free energy per spin as follows: 



3 The Methods 
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with K^'^^ representing the i-times transformed couphng constant K = J/ksT {ks - Boltzmann 
constant and T - temperature). The g in the above equation represents the contribution to the free 
energy (per spin) from degrees of freedom which have been decimated out in one RG transformation. 
For additional details of this method see for e.g. [T71 IHU]. 

As already stated the CSD will be used for the first time for quantum spin systems and we 
therefore describe it now in very detail. 

3.1 Configuration Selective Diagonalization 

As noted above exact diagonalization methods are limited to small systems because of an the expo- 
nential increase of the Hilbert space with the system size. In certain circumstances, however, the 
size of treatable systems may be increased significantly if only approximate energies and expectation 
values are required. Here we describe the adaptation of the CSD for the approximative calculation 
of low-energy properties of quantum spin systems in the context of the N=42 Sierpihski gasket. This 
approach has been developed using the methodology recently applied in the multi-configuration in- 
teraction methods fl^ El IHSj in quantum chemistry. In the following, we will describe this method 
and its implementation for the system at hand, but note that it is applicable to other systems as 
well. The division of the whole system into fragments which are exactly solvable within an exact 
Lanczos diagonalization is a necessary precondition for the application of the CSD, as discussed in 
detail below. The model in Fig |21 contains 42 sites, resulting in a Hilbert space of dimension 5.3 x 10^^ 



Figure 2: 42 Site approximant of the Sierpihski gasket. If the sites indicated by black circles are 
treated as a central fragment, the remainder of the system decomposes into three uncoupled thirteen- 
site satellite fragments. 

in the 5"*°* = subspace, which can be reduced by a factor of approximately four by discrete Abelian 
symmetries, such as spin- and real-space reflection. The use of the three-fold rotation incurs too 
much computational overhead to be of real value. Even so, it remains too large for exact diago- 
nalization schemes for low-lying eigenstates. The special structure of the cluster permits the use 
of approximate configuration-selective diagonalization methods. If we subdivide the cluster into a 
central fragment containing the three sites indicated with black dots, its remainder consists of three 
identical "satellite" fragments of 13 sites each. If there were no bonds between the central fragment 
and the satellites, the eigenstates of the overall system would be outer products of the eigenstates 
of the four fragments. In the presence of interactions, only four bonds couple the central fragment 
with each of the satellites, which may be assumed to perturb the spectrum of the satellites only 
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weakly. Because the satellite fragments are only weakly coupled to the central system, it is plau- 
sible to expand the eigenstate of the overall system in a basis of outer products of eigenstates of 
the satellite fragments. In the absence of any coupling to the central fragment, this approximation 
would be exact, in their presence one can hope that only a few configuration in this expansion will 
carry the overwhelming weight of the wavefunction. This observation suggests the application of the 
CSD (which has a long history in similar scenarios in quantum chemistry [211 ES|) to quantum spin 
systems. 

The basic idea of this approach is simple: Suppose an approximate wavefunction of the ground 
state is already known and this wavefunction has nonzero coefficients only for a small fraction of 
the configurations of the Hilbert space. We then estimate the weight of each remaining unselected 
configuration in second order perturbation theory. We keep only such configurations where the 
absolute value of the estimated coefficient surpasses a predefined threshold e and sum the perturbative 
energy contributions of the neglected configurations. Next we determine the eigenstate within the 
new selected Hilbert space with a direct diagonalization technique. The resulting state will be a better 
approximation of the desired eigenstate in the full Hilbert space. These two steps are alternated with 
decreasing selection threshold e and the energy (including the perturbative correction of the discarded 
configurations) and other expectation values are extrapolated to the limit e — > 0. For many systems 
this limit can be safely extrapolated with selected Hilbert spaces that contain only a small fraction of 
the possible configurations. The process is initiated with some simple trial wavefunction containing 
only the appropriate ground-state configurations of the segments. 

The algorithm thus consists of two distinct phases: in the expansion step new configurations are 
selected perturbatively and in the diagonalization step the lowest (or a few of the lowest) eigenvalues 
of the selected Hilbert space are determined. In the diagonalization step we iteratively improve a trial 
vector for the ground state of the selected Hilbert space using a preconditioned Davidson method 
jSH ES] • The time-consuming step of this iterative method is the computation of expectation values 
of the many body Hamiltonian H between trial states 

k 

where designate the configurations of the selected Hilbert space. The evaluation of such matrix 
elements is difficult, because at any given stage, the selected Hilbert space contains an essentially 
random subset of the possible configurations. 

The Hamilton operator H of the system can be written as 

H = H, + J2H, + Y.Hsc (4) 

s s 

where s enumerates the satellite fragments and c designates the central fragment. Hg and sum 
terms of H acting on a single fragment, while if^c couples the satellite s to the central fragment. 

To evaluate the expectation values we pre-diagonalize the 13-site satellite fragments in their 
respective Sl°^ spin-segments and compute the boundary-operators S*^ and 5"^ for the four boundary 
sites b of each fragment s in this basis. Similarly we compute the matrix representation of the 
corresponding operators for the central fragment. Each configuration |0fc)in the Hilbert space is 
labeled by a quadruplet of quantum numbers (rzc, ni, n2, ns), where is the index of an eigenstate 
of the corresponding fragment. The first two terms in the Hamiltonian Q are diagonal in this 
representation and easily evaluated. 

Nondiagonal terms, in contrast, are difficult to evaluate because of the sparsity of the selected 
configurations in the overall Hilbert space. In order to avoid costly lookup operations we have devel- 
oped a so-called residue based scheme for the evaluation of the matrix elements that we will describe 



6 



n3 



CI 



nlk(R,Cl,Sl) n2k(R,Cl,S2) 



C2 



C3 



C4 



C5 



C6 



Figure 3: A single entry R in the residue tree (top left) is designated by the quantum numbers of the 
fragment not changed by the operator in question (see text, {n2-,n^) in the example). The first level 
labels the quantum number of the central fragment that is necessary to completely specify the state. 
Attached to each such entry Cj is the list of indices of the configurations that have quantum numbers 
(cj, . . . ,n2,n3). Each element of such a list contains the quantum number rii of the configuration 
and its associated index of the coefficient of the associated configurations k{R, Ci,s). 



in detail in the following. Each coupling term Hgc is a sum of products of pairs of boundary operators 
described above. To efficiently evaluate the expectation values of this part of the Hamiltonian we 
use a residue based matrix element evaluation technique that was originally developed for selecting 
configurations in interaction methods fl^ 1ST] . For a particular configuration of the selected Hilbert 
space one individual term in Hcs changes the quantum numbers on the central fragment and on one 
of the satellite fragments (in the following we choose without loss of generality fragment 1). The 
quantum numbers on the other two fragments n2 , are the same on the right- and left-hand side of 
the configuration for nonzero matrix elements. All all nonzero contributions arising from this par- 
ticular term in the Hamiltonian connect two configurations which have the same index ^2, ?t-3. Since 
we wish to enumerate all nonzero matrix elements, we can use this set as a label for the associated 
matrix elements. In the following we will call this set of quantum numbers of both configurations 
the transition residue (R) mediating the matrix element. For each transition residue (^2, ''^3) we can 
enumerate the configurations that lead to nonzero matrix elements by a doubly nested residue tree. 

All nonzero matrix elements mediated by a particular transition residue R are uniquely enumer- 
ated by pairs of indices of the first level and pairs of indices of the second level of the residue tree R. 
There are no matrix elements spanning different residue trees. As an example the matrix element 
evaluation for Hsc = S^Si can thus be written as: 

inHcsm) = E E 

_R=(n2,n3) ci,C2in Tree R 

^ciC2 ( E E 

Slin list {R,ei) S2iii list (R,C2) 

^(i,Ci,Si)'5'siS2^(i,C2,S2)) (^) 
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In order to carry out the sums in the above equation, we construct a residue tree for each 
R = {^2,^^), which is illustrated schematically in FigOl The first level of the tree enumerates 
the allowed quantum numbers of the central fragment {ric = Ci, 02,03, . . .). The second level of the 
tree enumerates for each central fragment quantum number the selected configurations (q, Sj, n2, n^) 
and contains the value of the coefficient of the associated configuration: A{R = (^2, na), q, Sj, ). The 
sum in Eqn.(j3I) is then computed by picking all allowed pairs of branches Cj, Cj of the tree. For each 
pair of branches, each pair of entries generates a nonzero matrix element. This double sum can thus 
be performed without further lookup operations and evaluates nevertheless only matrix elements 
that connect selected configurations. Note that the innermost loops run over all selected indices on 
the satellites, i.e. encode 0(iVf ) operations for a fully selected Hilbert space, where A^^ ~ 1000 is the 
dimension of the Hilbert space on a satellite fragment. In order to build the tree just once for many 
applications of the Hamiltonian, it is more efficient to store the indices k{R, 02,82) rather than the 
coefficients in the tree. 

Using this procedure all matrix elements can be evaluated for arbitrarily complicated subspaces 
without lookup operations. We have implemented the residue tree by nested Adelson-Vleski-Landes 
(AVL) balanced binary trees which permit 0(log(N)) read/write operations per access. The numerical 
effort associated with building the residue tree is then proportional to the number of configurations. 
The number of matrix elements encoded by the tree, however is proportional to the expectation value 
of the square of the lengths of the inner subtrees, i.e. the inner sum in Eqn.(0). This sum scales with 
the number of configurations on a satellite fragment {0{Ng) for the fully selected case. As a result 
this matrix evaluation scheme is very efficient, in tests the expectation values were computed 200 
times faster than with a traditional hash-table based implementation. This increase in the numerical 
efficiency permits the treatment of much larger Hilbert spaces. 

Since the different residue trees are independent of one another, we can implement a relatively 
simple, scalable parallelization of the matrix element evaluation on a limited number of nodes, by 
distributing the residue trees across the nodes. In our implementation using 8 nodes of an SGI Power 
Challenge incurred a total loss of about 12 % of CPU time compared with a run on a single node. 
This loss of efficiency results from the fact that the number of matrix elements encoded by a given 
residue tree depends on the number of configurations containing its residue quantum numbers. 

In the expansion step for a given reference state l^') we diagonalize 



for each trial configuration and use its coefficient for the selection criterion. If the new configura- 
tions were mutually non-interacting this would generate the exact eigenstate of the system. The 
time consuming step in the diagonalization of equation ® is the computation of the off-diagonal 
coupling matrix element (0j|if|\E') between the reference state and the trial configuration, which can 
be accomplished using a similar residue based scheme as for the evaluation of the matrix elements. 
Initially we choose a small subset of configurations with the lowest energies which is diagonalized 
exactly. The lowest eigenstate of this Hilbert space is used as the reference state for the first iteration, 
for subsequent iterations the converged state of the previous iteration is used. 

For the present system, however, it is not feasible to even consider all possible trial configurations 
(N=O(10^^)) in the expansion loop. In each expansion step we therefore proceed as follows: first 
we order the configurations of the last state by the absolute weight of their coefficients. We then 
generate the interacting configurations for the most important configurations until the total weight 
of the configurations considered in the reference state exceeds about 90% of its norm. Interacting 
configurations of a given configuration are all those that have a nonzero matrix element with the 
Hilbert space of the selected configurations. These are typically only a small fraction (about 1-3%) 
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of the selected configurations. For these configurations we apply the selection criteria and gather the 
selected subset containing ni new configurations. We then generate the interacting configurations 
of the next important segment of the reference state until to total weight of the configurations of 
the reference state exceeds 96% of its norm and generate n2 new configurations from this set. We 
continue this process with geometrically decreasing fraction of the norm until the number of newly 
generated configurations Uk is less than 5% of the total number of newly generated configurations 
'^i + ^^-2 + ■ ■ ■ + ^fe-i- Since configurations with very little weight in the original wavefunction are 
unlikely to generate interacting configurations that will be selected we avoid to even build the full 
interacting space of the reference wavefunction. 

FiglUillustrates the convergence of the ground state energy with respect to the selection threshold 
for the system at hand. It demonstrates that for selection thresholds less then 1 x 10~^ the energy 
is extrapolated with an accuracy of one percent or better. 




Selection Threshold 



Figure 4: Convergence of the energy as a function of the selection threshold 

Expectation values, such as correlation functions can be easily computed by adapting the residue 
driven scheme. One constructs residue trees corresponding to the appropriate operators and evaluates 
the expectation values similar to the energy. 

4 Results 

4.1 Spin-spin correlations 

We turn now to the investigation of the ground state magnetic order on the quantum s=l/2 anti- 
ferromagnetic Sierpihski gasket. The careful investigation of the spin-spin correlations (SiSj) will 
provide a deeper insight into the ordering behavior. 

We begin with the ground state energy (being just the sum over all nearest-neighbor spin-spin 
correlations on the Sierpihski gasket). In the following table we provide the calculated value of the 
ground state energy for all finite lattices up to N=42. We provide the site and the bond average 
energies (This might help other groups to compare to our data.) In a quantum spin system with 
Neel-like semi-classical ordering the spin-spin correlation between the spins in one classical sublattice 
remains constant for large distances after a slight decay due to the quantum fluctuations. Therefore 
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Table 1: Number of spins N, site energy and bond energy Cf, for the Sierpinski gaskets up to N=42. 
(For N=42 the error of the CSD calculation is given in brackets.) 



N 


Heisenberg 




XY 




es 


eb 


es 


eb 


6 


-0.375 


-0.25 


-0.25 


-0.166666 


15 


-0.416125 


-0.231181 


-0.282024 


-0.156680 


42 


-0.439(5) 


-0.227(9) 


-0.298(8) 


-0.154(9) 



we first investigate the spin-spin correlations between the spins of one classical sublattice as defined 
in Sec El (for example all sites with circles • in Fig^. 

The geometrical distance, which is normally the measure on one- or two-dimensional lattices, is 
not easily transferable to fractal objects, such as the Sierpinski gasket. For the Sierpinski gasket 
we therefore use the Manhattan distance tm which counts the minimal number of steps required to 
connect one site to the other. Because of lack of translational symmetry in the Sierpinski gasket we 
find different spin-spin correlations between spins having the same Manhattan distance. In what 
follows we use a simple averaging procedure over all spin-spin correlations at a given Manhattan 
distance. 

In our earlier investigations we predicted a disordered ground state using the exact diagonalization 
data for N=15 among others. Below we will compare this data with the new data for N=42 to 
confirm our prediction. In FigElthe spin-spin correlation {SiSj)^j^^ for the Sierpinski gasket (SG) 
with N=42 and corresponding data for a two-dimensional square lattice (SL) with N=40 [2S1 12Z] 
(where magnetic long-range order is well known to exist) is shown. For the SL data we have chosen 
the shell distance rs (where shell-like circles numbered 1,2,3,. . . are drawn around a given site which 
then connect all neighbors at the given distance). The Sierpinski gasket shows a dramatic drop in 

0.25 I — > — > — > — > — > — > — > — > — I I — > — > — > — > — > — > — > — I 




Figure 5: The spin-spin correlation between spins on a classical sublattice {SiSj)^j^^ vs. distance. 
Left: A = 1 (Heisenberg model) - data for the Sierpinski gasket with N=42 (SG) vs. Manhattan 
distance vm and corresponding data for a two-dimensional square lattices with N=40 (SL) vs. the 
shell distance rs] Right: A = (XY model) - data for the Sierpinski gasket with N=42 for the x- 
and z- components of (SiSj). 

the spin-spin correlations and for large distances it even changes its sign. This is in marked contrast 
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to the two-dimensional square lattice where after a slight decrease due to the quantum fluctuations 
an almost constant behavior over distance is observed. This points to a complete loss of the classical 
Neel-like magnetic order in the quantum Heisenberg AFM on the Sierpihski gasket. Almost the same 
behavior can be seen for the XY model with A = 0. Here we have to distinguish between the different 
components of the spin-spin correlation due to the spin anisotropy: (S^Sj) = {SfSj) ^ {S-Sj). But 
as one can see the different components of the spin-spin correlations show similar behavior and the 
complete loss of any classical Neel-like magnetic order. 

Even though the classical magnetic order seems to be absent, some other type of long-range order 
in the pair correlations might prevail in the Sierpihski gasket. In order to check this conjecture we 
calculate the absolute spin-spin correlation | {SiSj)\ - j^j^ between all spins (and not just between spins 
in a classical sublattices) over vm (again using an averaging procedure as described above) and show 
a semilogarithmic plot of the data. From FiglHlwe deduce that only very short ranged correlations 




Figure 6: The average of all absolute spin-spin correlations \{SiSj)\ - j^j^ vs. Manhattan distance 

tm of the Sierpihski gasket. Left: A = 1 (Heisenberg model) - K^iSj)! for N=15 and N=42; Right: 
A = (XY model) - \ {SfSj) \ and \ {SfS]) \ for N=42. 

exist at all in the Sierpihski gasket and this behavior is independent of the spin anisotropy. For 
the XY model (right) we observe again a similar behavior for the different components of the spin- 
spin correlation. From the semilogarithmic plot we derive a correlation length ^ ^ 1, applying 
K^jSj)! ~ e"'-'"*^/^-'. The very small (in absolute values) upturn of the spin-spin correlations for the 
largest separation is presumably a boundary effect of the corner spins of the Sierpihski gasket. These 
corner spins have a different coordination number (only 2 bonds instead of 4) and they are mainly 
contributing to the value of this particular spin-spin correlation. 

As we observe only short range order in the Sierpihski gasket, we turn to the investigation of 
the local order which might exist in this lattice. We show in Fig|7|the spin-spin correlations on the 
Sierpihski gasket bonds only and chose (S^Sj) (Heisenberg case) and {SfSj) (XY case). The other 
components behave quite similar. 

We observe in the Heisenberg case (left) as well as in the XY case (right) a tendency to a plaquette 
formation at the corner spins. An example of such a plaquette is seen between the spins 39-40-42-41. 
Here the spin-spin correlations are very pronounced between the neighbor spins on the assumed 
plaquette (and reach about 80% of the true isolated plaquette value as we have checked) and they 
are apparently much smaller to the remaining lattice (where for an isolated plaquette we would 
observe 0). And even though there is an antiferromagnetic bond between the spins 40 and 41, the 
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Figure 7: The nearest neighbor spin-spin correlation on the Sierpihski gasket with N=42. Left: 
A = 1 (Heisenberg model) - {S^sf), Right: A = (XY model) - {SfSp (The dotted lines stand for 
a ferromagnetic correlation). 

resulting spin-spin correlation between this two spins is ferromagnetic which points to a beginning 
dimer formation between them. The other dimer correlation (S'39S'42) (value not shown in Fig|7j) is 
0.234 (Heisenberg case) and 0.230 (XY case) and therefore very close to a true dimer correlation of 
1/4. One could argue that another plaquette may form inside the lattice (one example might be 
spins 30-31-34-33). We have found a similar behavior already for the N=15 Sierpihski gasket. The 
building of plaquettes in a lattice with strong frustration seems quite interesting given the fact that 
in the strong frustrated region of the Ji — J2 square lattice model one might find a similar behavior 
(although this behavior is still under controversial discussion) jSEl IHZl IHB] • 

4.2 Spin gap 

In this section we investigate the spin gap F, defined as 

T = E{S^i^ + 1) - E{S^iJ. (7) 

The absence of Neel-like magnetic long-range order is as a rule accompanied by a finite spin gap and 
can therefore be used as another criterion for a disordered ground state jSHl- In FiglSlwe present 
F for the Heisenberg and XY case. The data for N=6 and 15 has been presented before but with 
the CSD we are able to calculate the spin gap for N=42 as well. The data has to be analyzed with 
particular care, because we expect finite size effects still to be present. If we just make a simple linear 
fit through the data we roughly get a value F ^ 0.2. One might also argue that N=6 is too small 
to be taken into account for this consideration and therefore we did another fit with only N=15 and 
42. Still we see that the spin gap F remains finite, but is reduced to F ^ 0.1. A similar conclusion 
holds for the XY model, although the spin gap appears to be smaller (this behavior has been found 
in other investigations too [ID])- Though our finite-size extrapolation must be taken with particular 
care we see further arguments in favor of a finite spin gap and a ground state without magnetic 
long-range order. 

We mention that the AFM on the Sierpihski gasket belongs to the class of frustrated spin systems 
(like the kagome or the checkerboard lattices) having exactly known localized magnon eigenstates 
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Figure 8: The spin gap vs. inverse system size for the Sierpihski gasket with N=6,15 and 42. Left: 
A = 1 (Heisenberg model); Right: A = (XY model). 

leading to a macroscopic magnetization jump to saturation j^T] |32]- These locahzed magnons can 
live e.g. on 'hexagons' (large triangles) inside satelhte fragment (e.g. the sites 4,5,6,12,11,8 in Fig|7j). 
For N=42 the magnetization jumps at the saturation field hgat = 3J from m = 17/21 to saturation 
m = 1. Even for N — oo the height of the jump remains finite since the number of localized magnons 
occupying the lattice growths with N. 

4.3 Low temperature thermodynamics 

It is known that the low-temperature specific heat is closely related to the low-lying excitations of a 
system jH 021 Ell IS] ■ As we have argued already in Sec 14. 21 the low-lying excitations in turn might 
show a spin gap behavior and therefore point to a disordered ground state. We and others have 
seen in previous investigations a close connection between an additional low-temperature peak in 
the specific heat and a finite spin gap |18t |S1 HEl El • Following this argumentation we will analyze 
the specific heat Cy of the system especially in its low-temperature region. We show in Figj^l c„ for 
the Heisenberg and XY Sierpihski gasket calculated with CD and QDT. We note that the c„ results 
for the Heisenberg case are identical to those of Ref. and shown for comparison purposes. We 
observe in both cases additional low-temperature peaks which relate to two different energy scales 
relevant in the system. The first energy scale is connected to the typical broad peak, whereas the 
second one is connected to the low-temperature peak and its value is connected to the finite spin 
gap Foo- This behavior of the specific heat is typical for all antiferromagnetic systems on "corner 
sharing triangles" lattices {kagome Sierpinski ^71, squagome [STF). In fact, in all those systems 
the basic unit leading to this behavior is a A-chain jJUl which shares spins with other A-chains. 
{kagome - 12 spin chain, Sierpinski - 6 spin chain, squagome - 8 spin chain). The RG transformation 
used here takes into account the excitation spectrum of 6 spin chain (Fig. [H right) and by using 
Eqn.Q one gets a larger low-temperature peak in comparison to exact diagonalization data. 

In both cases an additional low-temperature peak in the specific heat constitutes an additional 
argument for a finite spin gap and therefore for a disordered ground state. The small additional 
peak in the exact diagonalization data at even lower temperatures has been attributed to a finite 
size effect P^ . 
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Figure 9: The low temperature specific heat: full diagonalization data for the N=15 system and 
quantum decimation data. Left: A = 1 (Heisenberg model); Right: A = (XY model). 



5 Summary 

We have carried out a numerical investigation of a s= 1/2 quantum antiferromagnet on the Sierpihski 
gasket with two types of spin exchange, Heisenberg and XY. We have used the exact diagonalization, 
a newly implemented configuration selective diagonalization approach and a quantum decimation 
technique to calculate the spin-spin correlations in the ground state, the spin gap and the low- 
temperature specific heat. 

The current investigation complements and verifies previous work done for the Sierpihski gas- 
ket |15| ITHl [T71 ITHl [T^ 1^ . Main progress in the investigation of this fractal many-body system 
results from the successful introduction of the configuration selective diagonalization (CSD). This 
new method permits the calculation of the ground state wavefunction (and therefore all correlation 
functions) and excited states as well for larger finite quantum spin system. Using CSD we calculated 
the spin-spin correlations and the spin gap for N=42. We note that the method can be applied to 
other frustrated low-dimensional quantum spin systems with just moderate changes. 

The reported data suggest that the ground state of the Sierpihski gasket remains disordered for 
Heisenberg and XY spin exchange. It seems that the interplay of quantum fluctuations and low 
dimension prevents any kind of magnetic long range order in this system. From the available data 
we derive a magnetic correlation length ^ ~ 1. The nearest-neighbor spin-spin correlations show a 
tendency to plaquette formation. It will be interesting to find out whether or not this behavior is 
related to similar findings in the strongly frustrated region of the Ji — J2 Heisenberg antiferromagnet. 
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